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The dynamical behavior of matter wave solitons of two-component Bose-Einstein condensates 
(BEC) in combined linear and nonlinear optical lattices (OLs) is investigated. In particular, the 
dependence of the frequency of the oscillating dynamics resulting from initially slightly displaced 
components is investigated both analytically, by means of a variational effective potential approach 
for the reduced collective coordinate dynamics of the soliton, and numerically, by direct integrations 
of the mean field equations of the BEC mixture. We show that for small initial displacements 
binary solitons can be viewed as point masses connected by elastic springs of strengths related to 
the amplitude of the OL and to the intra and inter-species interactions. Analytical expressions of 
symmetric and anti-symmetric mode frequencies, are derived and occurrence of beatings phenomena 
in the displaced dynamics is predicted. These expressions are shown to give a very good estimation 
of the oscillation frequencies for different values of the intra-species interatomic scattering length, 
as confirmed by direct numerical integrations of the mean field Gross-Pitaevskii equations (GPE) of 
the mixture. The possibility to use displaced dynamics for indirect measurements of BEC mixture 
characteristics such as number of atoms and interatomic interactions is also suggested. 

PACS numbers: 03.75.-b,67.85.Hj, 05.45. Yv 



I. INTRODUCTION 



Binary mixtures of Bose-Einstein condensates (BECs) 
are presently attracting a great deal of interest in connec- 
tion with a series of interesting phenomena such as the 
formation of segregate domains [l[ , polarized states , 
spin textures [^^topological excitations [Jl, novel Joseph- 
son oscillations [a, Q Rabi Josephson oscillations Q , four 
wave mixing [8|, etc. Moreover, multi-component BECs 
have been shown to support nonlinear waves of novel t ype 
such as symbiotic solitons [9|, domain- wall solitons [lO[ 
and combinations of dark-dark [ll| and bright-dark soli- 
tons [l2, [13 1 the last one leading to long lived oscilla- 
tions which were experimentally observed in |14{ . The 
possibility to trap binary mixtures in optical lattices 
(OLs), experimentally demonstrated in [l5|, has added 
further interest to the field. In particular, the interplay 
between the nonlinearity induced by the interatomic in- 
teractions and the strength of the OL has been shown to 
lead to interesting phenomena such as Landau- Zener tun- 
neli ng [ l6j. transitions from superfluids to Mott insula- 
tors [l7|. Moreover, existence of nonlinear periodic waves 
on nonzero backgrounds [l^l, gap solitons [l^l, mixed- 
symmetry modes and breathers both in continuous and 
discrete (arrays) mixtures [20|. 

Besides usual (e.g. linear) OLs, it is also possible to in- 
troduce a periodic structure in the system by modulating 
the scattering lengths in space by means of the Feshbach 
resonance technique [21|. This allows to create what is 
known as a nonlinear optical lattice (NOL). BEC mix- 
tures in NOLs have been recently considered in connec- 
tion with quantum simulation of novel Hubbard models 
[22l | and interesting phenomena such as sonic analogues 



of black holes l23| and control of soliton creation [2J]. 
A possibility of observing delocalizing transition even in 
one-dimensional BECs loaded in OLs due to the presence 
of the NOL has been also suggested [25| . For a fresh re- 
view on BECs in nonlinear optical lattices we refer the 
article in [2g|. In all these studies, however, the effects 
of a combined linear and nonlinear optical lattice on the 
soliton dynamics and the link between dynamical behav- 
iors and interactions, have been scarcely investigated. 

The aim of the present paper is to study the mean 
field dynamics of initially displaced soliton components 
of binary BEC mixtures in the presence of a combined 
linear and nonlinear OL. In particular, the dependence of 
the frequency of the resulting oscillating dynamics on the 
inter-species interaction and on the number of atoms is 
investigated. This is done both analytically, by means of 
a variational effective potential for the displaced dynam- 
ics, and numerically by direct integrations of the mean 
field equations of the BEC mixture. We show that in 
the limit of small initial displacements, the effective po- 
tential leads to a mechanical interpretation of a binary 
soliton motion in terms of two point masses connected 
by elastic springs of strengths related to OL's amplitude 
and to the intra and inter-species interactions. The dis- 
placed dynamics, being the same as the one of coupled 
harmonic oscillators, can be decomposed in term of a nor- 
mal mode analysis from which analytical expressions of 
the symmetric and anti-symmetric mode frequencies, are 
explicitly derived. These expressions are shown to give 
a very good estimation of the oscillation frequencies for 
different values of the intra-species interatomic scatter- 
ing length, as confirmed by direct numerical integrations 
of the mean field Gross-Pitaevskii equations (GPE) of 



the mixture. The occurrence of beating phenomena for 
unequal and for equal numbers of atoms in the mixture 
for small interspecies interactions, is also discussed. The 
stabilities of stationary and oscillating dynamics are in- 
vestigated by Vakhitov-Kolokolov (VK) criterion [271 and 
by numerical simulations, respectively. These results sug- 
gest the possibility to use dynamical behaviors of suitably 
prepared initial multi-component BEC solitons as a tool 
for extracting information about physical characteristics 
of BEC mixtures such as interatomic interactions and 
species populations. 

The paper is organized as follows. In Sec. II we in- 
troduce the mean field model equations describing BEC 
mixtures in combined linear and nonlinear optical lat- 
tices and derive a variational effective potential formula- 
tion for the matter waves soliton dynamics. In section III 
we consider the displaced binary soliton dynamics in the 
framework of a coupled harmonic oscillator model which 
is valid in the limit of small initial displacements. Ana- 
lytical expressions for the symmetric and antisymmetric 
mode frequencies are explicitly derived. In Sec. IV re- 
sults of displaced soliton dynamics obtained by direct nu- 
merical integrations of the CPE are compared with the 
analytical predictions. The stability of stationary two 
component solitons and their slightly displaced dynam- 
ics are also investigated. Finally, in Sec. V the main 
results of the paper are briefly summarized. 



II. MODEL EQUATION AND VARIATIONAL 
ANALYSIS 



We consider as a mean field model for a mixture of two 
homonuclear condensates [28| in an external trapping po- 
tential, the following system of coupled Gross-Pitaevskii 
equations 



dt 



2m dx^ 



(1) 



where ipj (j ~ 1,2) denote the wave function of the binary 
mixture and Vext {x) the external potential resulting from 
harmonic and optical lattice confinement, in the following 
taken of the form 



Vext{x) = -mwlx'^ + Vl cos{2kLx). 



(2) 



Here uJx and uj± are the longitudinal and transverse trap- 
ping frequencies of the harmonic confinement, as and 
tts are the intra- and inter-species scattering lengths, 
Vl and fci are respectively strength and wave number 
of the optical lattice. Since the longitudinal harmonic 
confinement introduces only slight modifications to the 
longitudinal periodic potential (in experimental settings 
LOx is of the order of a few Hz), it will be ignored in the 
following [29|. Introducing dimensionless variables: 
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Eqs. [T]can be written in the form 
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Vl 

"FT 
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rcscaled intra- and inter-species interaction strengths. In 
this Eqs. the order parameter ipj is normalized to the to- 
tal number of atoms such that /_^ (l'0iP + l'/'2p) ds = 
^1 +-^2; where Nj, j = 1,2 are the separately conserved 
numbers of atoms in each component. In the following we 
fix fci = 2 and assume a dependence of the intra-species 
interaction of the form 



(0) , (1) /r, N 

5ii = .9ii + .9ii cos(2s) 



(5) 



with the spatial modulation part denoting a NOL of 
strength Pii . In an experimental context such a spa- 
tial modulation could be produced by optically induced 
Fcshbach resonances [30| e.g. by a laser field tuned near 
a photo association transition. Virtual radiative tran- 
sitions of a pair of interacting atoms to this level can 
then change the value and even reverse the sign of the 
scattering length. It can be shown that a modulation 
of the laser field intensity of the form / = /oCos^(kx) 
reflects in a modulation of the scattering length of the 



.(1) 



form aV\x) = ^^^[l + al / {6 -f /)] 



where a^g is the 
intra-species scattering length in the absence of light, & 
is the frequency detuning of the light from the resonance, 
and a is a constant factor [30l . l31j . For weak intensities 
Iq <C \5\ the real part of the scattering length can be 
then approximated as Os = a^g -l-a^j^ cos^(kx) which is 
essentially the same form assumed in Eq. ([5]). 

Note that in the absence of the OLs and with gi2 = 
0; Eq. ^ decouple into two nonlinear Schrodinger 
equations which admit, for attractive intra-species in- 
teractions, exact bright soliton solutions with typical 
Gaussian-like function shape. With the view to solve 
Eqs. in ^ within a variational approach, we then adopt 
for the coupled soliton wavefunction the following ansatz 

(g — 5„ ■)2 

V'j(s,t) = yljcxp[ ^ 2 +»(^Oj(s-soj) + ^j)], 



2ai 



J = 1,2 



(6) 



with parameters Aj 



SQj, $j, denoting amplitude. 



width, center of mass and phase of the soliton, respec- 
tively, taken in the following as time-dependent param- 
eters. Note that the wavefunction is normalized to the 



total number of atoms Nj so that Aj 



N, 



^/ttuj ' 



The effective Lagrangian for the system is written as 
(£) = J_ Lds with the Lagrangian density L given by 
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From the Ritz optimization conditions [32|, we have 
^ = 0, ^ = 0,^ = and '^ ^0. The first 
optimization condition 
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0, 



(8) 



in conjunction with the normaUzation condition of ipj 
impUes that ^/nujA'^ = Nj is a constant. This constrain 
when used in the relations obtained from the other opti- 
mization conchtions give 
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and 
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In order to derive an exphcit formula for the effective 
interacting potential of the coupled solitons, we consider 
that the condensates are symmetrically placed with re- 
spect to the a OL minimum i.e, soj = iso/2- In this 
case, Eqs. in (fTTj) can be combined to give 



sin(so) 



so - 2Va (^e-'^i + 
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FIG. 1: Effective potential versus separation for Vo = —0.5. 
Top panel gives Vcft with A'^i = 1.0 and A'^2 ~ 0.5 for dif- 
ferent values of gi2, namely, —0.2 (solid), —0.4 (dotted) and 
— 0.6(dashed). The middle panel gives Vcff with A''i = 1.0 
and A'^2 = 1.0 for different values of 312, namely, —0.2 (solid), 
—0.4 (dotted) and — 0.6(dashed). The bottom panel shows 
14ff with 312 = —0.5 for different values of TV = A^i = N2, 
namely, iV = 0.4 (solid), 0.8 (dotted) and 1.2 (dashed). 



as the evolution equation for the separation sq between 
center of the solitons. Notice that Eq. (fT2)) is the same 
as the dynamics of a Newtonian particle in the effective 
potential 



l/eff(so) = [2^0 (e~"' + e-"^') cos(.so) 
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(13) 



Also notice that this potential has the absolute mini- 
mum in the origin and that for small values of sq around 
the minimum of the potential can be approximated as a 
harmonic potential. In such approximations, the small 
oscillation frequency of displaced solitons dynamics can 



be written as 
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Moreover, one can show that the vanishing condition 
of -tM- gives the chemical potential, /i, of stationary com- 
ponents as: 
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FIG. 2: Mechanical model of displaced binary soliton dynam- 
ics in terms of harmonic oscillators of elastic constant k cou- 
pled by a spring of elastic constant K12 . 
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(in writing Eq. (jisp we have used $j = —fJ-jT and soj = 
in Eq. ([S])). This expression can be used (sec below) to 
study the stability of stationary two component solitons 
through the Vakhitov-Kolokolov (VK) criterion. From 
Eq. (fT3|) we sec that the effective potential for the cou- 
pled solitons dynamics is highly anharmonic and consists 
of three terms: the first two arise from the linear and 
nonlinear optical lattices while the third one comes from 
the mutual interaction between the solitons. The mutual 
interaction term depends both on the number of atoms 
in the condensates and on the strength of the interac- 
tions. This part of the potential will therefore change 
sensitively with the variation of N and 512. 

In Fig. [T] we show the effective potential as a func- 
tion of So for two attractively interacting solitons and 
different values of —gi2 (left and middle panels) and 
N ^ Ni ~ A''2 (right panel). More specifically, left panel 
give Veff with A^i = 1 and N2 = 0.5 while middle panel 
shows Veff with Ni = N2 = 1 for different values of —512. 
Note that the inter-species interaction is effective mainly 
for BEC components with a significant spatial overlap- 
ping e.g. when they arc very close to each other. In this 
situation an oscillatory dynamics of the BEC components 
around their common center of mass can be started by 
slightly displacing them from the equilibrium position 
corresponding to the fundamental minimum of the effec- 
tive potential in Fig. [T] Also note that for an attractive 
inter-species interaction the absolute minimum of the ef- 
fective potential becomes deeper and deeper as iV|9i2| is 
increased. Therefore, the reduced equation of motion in 
(|T4l) implies that the solitons oscillation with respect to 
each other if they arc placed very close to the effective 
potential minimum at sg = 0. For repulsive inter-species 
interactions, however, the effective potential will have the 
shape of a barrier (rather than a potential well) with a 
maximum (rather than a minimum) at the origin. In this 
case, the soliton components move away from each other 
keeping their shapes unchanged |33| . 



III. NORMAL MODE ANALYSIS OF 
DISPLACED BINARY SOLITON DYNAMICS 

It is useful to gain some modeling insight of the dis- 
placed binary soliton dynamics in the limit of small dis- 
placements soi <C 1 and S02 ^ 1- In this case Eq. (fTTj) 
reduces to 



so, - iwoe-^^ + ^glk^^l^ 
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27ra,- 
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Let us concentrate for simplicity on binary solitons with 
equal number of atoms and equal widths e.g. A'^i = N2 = 
N and ai = 02 = a. Introducing parameters 
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we can rewrite Eq. (|16p in the form 
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(17) 
(18) 

(19) 

(20) 
(21) 



which are the same as the equation of motion of two cou- 
pled identical harmonic oscillators of mass M and elastic 
spring K connected by a spring of elastic constant K12 (see 
Fig. [3]). In the absence of inter-species interaction, (as 
it is the case, for example, when the interspecies scatter- 
ing length is detuned to zero by means of a Feshbach 
resonance) In the presence of interspecies interactions 
the above equations are readily decoupled in the nor- 
mal mode coordinates: fi = sqi — 502, C2 = soi + 502, 
this giving M^i = —uf^i, i = 1,2, with characteristic 
frequencies 



Wi = ± 



2ki2 



M 



W2 = ±4 / 



(22) 



and explicit normal mode solutions 



1,2. 



(23) 



The most general solution of the displaced soliton dy- 
namics in the coupled harmonic oscillator approximation 
follows from Eq. ([25]) as 



1 



1 



soi(i) = 7T(6W+6(i)), so2{t) = -iUt)~Ut))- (24) 



From these equations we see that the solution ^i asso- 



ciated to the frequency oji 



UJn 



corresponds to an 



asymmetric (out of phase) oscillation of the displaced two 
component soliton, while the solution ^2 corresponds to a 
symmetric (in phase) motion of frequency UJ2 = uJsym in 
which the coupling spring remains unstretched. Notice 



that I 



^asym 



is the same as the expression of the frequency 



derived in Eq. ([Ti)) . Also note that in analogy with opti- 
cal and acustical vibrations of molecules, this frequency, 
for attractive inter- and intra-species interactions, is al- 
ways higher than the frequency ujsym of the symmetric 
mode, e.g. uJasym/^sym > 1, with the equality hold- 
ing in the case 512 = 0. More explicitly, the following 
dependence for the frequency ratio of asymmetric and 
symmetric modes on parameters of the binary mixture, 
is derived: 
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Ngue'' /'^/a 
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(25) 



Note that in the weak coupling limit \gi2\ ^ 1 the dis- 
placed dynamics will display typical beating phenom- 
ena with a high frequency component oscillating in- 
side a slowly varying envelope, with beating frequencies 



^beat — ^asym ^sym^ 



plus order combinations. 



For the general case TVi 7^ N2 (equivalently, oi 7^ 02), 
the dependence of characteristic frequencies of the oscil- 
lators on parameters can be derived in similar manner. In 
the next section we shall compare the above predictions 
for the soliton displaced dynamics with direct numerical 
integrations of the coupled GPE in (|1|) . 



IV. DYNAMICS AND STABILITY OF 

DISPLACED BINARY SOLITONS: NUMERICAL 

RESULTS 

In the following coupled GPE numerical investigations 
we assume that solitons prepared in such a manner that 
their relative coordinates are located at small distances 
from the minimum of Veff in Fig. [TJ We remark that 
initial small displacements of the two components of 
the mixture could be experimentally induced by a rapid 
change of the inter-species scattering length from nega- 
tive to positive and then to negative again, by means of 
the Feshbach resonance technique. The inversion of the 




FIG. 3: Oscillation of coupled BEG components for Vo = 
—0.5, g{-i^' = —1, g{^ = —0.5 and sq — 0.2. Here, number of 
panels is counted from the top. First panel shows motion of 
soliton profiles for Ni — N2 = 1.0 and (712 = —0.2. Second 
panel. Same as in first panel but for Ni — N2 = 1.4 and 
gi2 = —0.5. Third panel. Motion of the centers of soliton 
components for unequal number of atoms A''i = 1, N2 = 0.5 
anf for gi2 — —0.2. The bottom panel show beatings aris- 
ing from the superposition of the oscillatory components dis- 
played in the 3rd panel. In all panels curves with big circles 
give results for PDEs in Q while dashed curves represent 
results for ODEs in (fTTj). 



sign of the interaction for a small fraction of time can 
be achieved with a properly designed time-dependent ex- 
ternal magnetic field. The component solitons will move 
in opposite directions during the short repulsive inter- 
species interaction time, and will become slightly sepa- 
rated (separation can be made small by properly reduc- 
ing the repulsive time). Taking into account detectable 
length scales of real experiments |34{ . we use in most of 
the calculations sq = 0.2. although larger initial displace- 
ments (sq ~ 1) will also be used for anharmonic effects. 
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FIG. 4: Beating dynamics of displaced binary BEG solitons 
with an equal number of atoms Ni = N2 = 1, for (712 = 
— 0.2, soi = 0.1, S02 = 1.3 (left panel) and 312 = — 0.5, soi = 
0.1, S02 = — 1.3 (right panel). Other parameters are fixed as 

.,(0) _ 1 „(i) 



ffir = -l>9li' = -0.5,Vo = -0.5. 



In Fig. [3] we show typical dynamics of displaced bi- 
nary solitons arising from a symmetric initial displace- 
ments with respect to the effective potential minimum. 
The top two panels refers to the case of unequal numbers 
of atoms. We see that in this case the soliton compo- 
nents oscillate with the same frequency which depends 
on inter-species interaction strength and on number of 
atoms in the condensates (compare 1st and 2nd panel). 
For Ni y^ N2, however, the oscillation frequencies of each 
component become unequal (3rd panel) with appearance 
of well-known beating phenomenon. The general solu- 
tion in Eq. ([M|) shows that the beating phenomena 
is expected also for equal number of atoms and small 
inter-species interactions if the motion is started with a 
generic initial displacement |soi| 7^ |so2|- This is exactly 
what the PDE calculations in Fig. 2] show for the case 
-^1 = -^2 = 1 and I Soil 7^ |so2| , in agreement with our 
normal mode analysis. 

The dependence of oscillation frequency (v) on gi2 is 
depicted in Fig. [5] for the case N1/N2 = 1. In particu- 
lar, the top left panel of this figure shows v vs —.912 for 
closely spaced solitons {sqj ~ ±0.1) while the top right 
panel gives the dependence on gi2 for sqj — ±0.5. No- 
tice that in both cases one can estimate values of v very 
close to the exact numerical results (dotted curve) from 
our simple analytical calculation. It may be pointed out 
that the variational approach gives little higher values of 
h' than that of the numerical integration for sqj = ±0.1( 
top panel) while we observed the opposite in the bottom 
panel. Also note that for overlapped condensates results 
obtained from analytical formula in P^ and ODEs in 
p2p are same. However, with the increase of separation 
deviation between these two results becomes appreciable. 
The deviation of analytical curve (solid) with ODE calcu- 
lations (dashed) implies that the effects of aharmonicity 
has been captured. In addition, the curves in both the 
panels clearly indicate that the values of ly increase as 
gi2 increases. This correlates with the fact that the cor- 
responding effective potentials become more deep and ac- 
quire larger curvatures at the origin as these parameters 
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FIG. 5: Top left panel. Frequency of the oscillation of the 
BEG components vs the inter-species interaction strength for 
fixed number of atoms A^i = A''2 = 1 and so = 0.2. Top 
right panel: Same as that in the left panel but for a larger 
initial displacement so = 1.0. Middle left panel. Frequency 
of the oscillation of the BEG components vs the inter-species 
interaction strength for the case of unequal number of atoms 
TVi = 1 and N2 = 0.5. and so = 0.2. Middle right panel: 
Same as that in corresponding left panel but for so = 1.0. 
Bottom left panel. Frequency of the oscillation of BEG com- 
ponents versus the number of atoms for fixed inter-species 
interaction strength gi2 — —0.5 and so = 0.2. Bottom right 
panel. Same as that in the left panel but so = 1.0. In all 
the panels the continuous curve represents the analytical ex- 
pression in equation (jldp . the dashed curve stands for result 
obtained from ODE in (|12p and open circles denote numerical 
GPE calculations. In each case, other parameters are fixed 
as: Vo = -0.5, g["^^ = -1 and g[\^ = -0.5 



are increased, clearly leading to higher frequency values. 
In the middle panel of Fig. [5] we portrayed a plot 
similar to that in the top panel but with TVi 7^ N2 and 
observed same behavior of the frequency curve. However, 
the observed frequency in this case is little smaller than 
the previous one. This might be associated with the fact 
that effective inter-species interaction is relatively larger 
Ni = N2 that Ni ^ N2 (Fig. [T]). In order to find de- 
pendence of v on number of atoms in the condensates. 
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FIG. 6: Relative frequency (fr = - 
and symmetric modes vs number of atoms N (left panel) and 
inter-species interaction gi2 (right panel). In both panels the 
solid curve refers to Eq.((25} while the open dots refer to GPEs 
numerical integrations. Other parameters are fixed as Vo = 
-0.5,ffl°' = -l,ff(i' = -0.5,so = 0.2. 
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FIG. 7: Time evolution of displaced component binary 
soliton densities. In each contour plots we have taken 
slightly perturbed initial conditions for parameter values Vo ~ 
—0.5, gjj^ = —1, (?jj = —0.5. Left and middle panels show 
density evolution for A^i = A'^2 = 1 and gi2 — —0.2 (left panel) 
and for Ni — N2 — 0.8 and 312 — —0.5. Right panel refers 
to the case of unequal number of atoms A''i = 1, A'^2 = 0.5 for 
512 = -0.2. 



we consider only the case N = Ni — N2 and plot in 
the bottom panels of Fig. [5] variation of i' with A^ for 
gi2 = —0.5. Particularly, the left bottom panel shows 
ly vs N for sq ~ 0.2 while the right bottom panel gives 
a similar plot but sq = l-O- As in the previous case, 
the frequency of oscillations here also increases with the 
increase of number of atoms due to the increase of effec- 
tive intcr-spccies interaction. Notice that the analytical 
(solid) and/or ODE (dashed) calculations here also gives 
a very good estimate of v of the exact numerical calcu- 
lation (dots with circles). The expected discrepancy of 
analytical curve from the ODE curve due to anharmonic 
effect is prominent for sq = 1-0 (bottom right panel). In 
Sec. Ill we showed that the frequency of the asymmetric 
mode is always greater than that of the symmetric one 
with a ratio Vr that depends on both 1712 and A^. In Fig. 
|6]we compare the dependence of Vr on 512 and on N as 
obtained from Eq. ([25|) with the one obtained from PDF 
calculations. We see that in both cases a relatively good 
agreement is found. It is also clear, both from analytical 
and numerical results, that Vr increases with the increase 
of either inter-species interaction and number of atoms, 
and that the analytical results are slightly overestimating 
this growth. 

The stability of the oscillatory soliton dynamics has 
been checked with slightly perturbed initial BEG profiles 
for given values of parameters and then allowed them 
to evolve according the GPE copupled equations. Den- 
sity plots for the evolution of soliton profiles for different 
values of 512 and number of atoms is displayed in Fig. 
[71 This figure clearly indicates that during time evolu- 
tion soliton profiles remain undistortcd. To confirm the 
stabile evolution of the density profiles we slightly vary 
initial conditions and noticed that they still evolve uni- 
formly with time. The stability of soliton profile can also 
be examined from the phase plot of coupled ordinary dif- 
ferential equations in (jlip . We have verified that in each 
case considered by us the phase plot exhibits stable focus. 



CONCLUSION 



In this paper we have studied the dynamics of matter 
wave solitons of two-component Bose-Einstein conden- 
sates in combined linear and nonlinear optical lattices. 

In particular, we have investigated the dependence 
of the oscillating dynamics resulting from two initially 
displaced BEG soliton components on the inter-species 
interaction and on the number of atoms. We showed 
that for small initial displacements binary solitons can 
be viewed as point masses connected by elastic springs 
of strengths related to the amplitude of the OL and to 
the intra and inter-species interactions. The displaced 
dynamics in can be decomposed in term of normal 
mode analysis from which analytical expressions of 
the symmetric and anti-symmetric mode frequencies, 
are derived. The occurrence of beating phenomena 
both for unequal and for equal numbers of atoms for 
small interspecies interactions, was also predicted. 
The stability of the oscillating dynamics has been also 
investigated by direct numerical GPE integrations. The 
predictions of the effective potential approach are found 
to be in quantitative agreement with direct numerical 
simulations. These results suggest the possibility to 
use dynamical behaviors of suitably prepared initial 
multi-component BEG solitons as a tool for extracting 
information about physical characteristics of BEG 
mixtures such as interatomic interactions and species 
populations. In this respect, we remark that in contrast 
to intra-species interactions, direct measurements of 
the inter-species scattering lengths are more difficult to 
access. The possibility to measure interspecies scatter- 
ing lengths through dynamical behaviors of displaced 
BEG components represents therefore an interesting 
possibility which could be tested in real experiments. 
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